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Abstract. We provide a new representation of a refinable shift invariant space 
with a compactly supported generator, in terms of functions with a special 
property of homogeneity. In particular these functions include all the homo- 
geneous polynomials that are reproducible by the generator, what links this 
representation to the accuracy of the space. We completely characterize the 
class of homogeneous functions in the space and show that they can repro- 
duce the generator. As a result we conclude that the homogeneous functions 
can be constructed from the vectors associated to the spectrum of the scale 
matrix (a finite square matrix with entries from the mask of the generator). 
Furthermore, we prove that the kernel of the transition operator has the same 
dimension than the kernel of this finite matrix. This relation provides an 
easy test for the linear independence of the integer translates of the genera- 
tor. This could be potentially useful in applications to approximation theory, 
wavelet theory and sampling. 



1. Introduction 
A function ip : R — > C is called refinable if it satisfies the equation: 

N 

(1.1) ip(x) = ^c k ip(2x - k), 

for some complex scalars Co, ...,cn. The scalars c k are the mask of the refinable 
function. We consider the case in which tp is compactly supported. Define the Shift 
Invariant Space (SIS) generated by tp as: 

S(tp) = {f : R — ► C : f{x) = £ y k tp{x + k), y k € C}. 

A refinable SIS is a SIS with a refinable generator. Refinable SIS and refinable 
generators have been studied extensively, since they are very important in Approx- 
imation Theory and Wavelet Theory. 

Many properties of ip can be obtained imposing conditions on the mask. One 
fundamental question is when the space S(ip) contains polynomials and of which 
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degree. The accuracy of tp is the maximum integer n such that all the polynomials 
of degree less or equal than n — 1 are contained in S{tp). 

The accuracy is related to the approximation order of S(tp), ( Jia95 , dB90 
and references therein), and with the zero moments and the smoothness of the 
associated wavelet when tp generates a Multiresolution Analysis Mey92| . There 
are many well known equivalent conditions for accuracy. The one that interests us 
here is the following, [Dau92| . CHM98 : 

Proposition 1. Let ip be a compactly supported function satisfying Then tp 

has accuracy n if and only if {1, 2 _1 , 2~(™ -1 ))} are eigenvalues of the (N + 1) x 
(N + l) matrixT defined by T = {c2i-j}i 1 j=o,...,N ! (the scale matrix,) and there exist 

polynomials po, ...,p„_i with degree{pi) = i such that the vector Vi — {pi(k)}k=o n 

is a left eigenvector of T corresponding to the eigenvalue 2~ l . 

Here, and always throughout the paper, we assume Ct = 0, if t ^ 0, . . . , N. One 
interesting property is that if tp has accuracy n, then for s = 0, 1, n — 1 it is 
true that x s = ^2 kG ^ J Ps(k)tp(x — k), where p s is the polynomial that provides the 
eigenvector for the eigenvalue 2~ s . So, if the polynomial x s is in S(<p), then the 
left eigenvector of T, corresponding to the eigenvalue 2~ s , provides the coefficients 
needed to write x B as a linear combination of the translates of tp. 

A local basis for S(tp) is a set of functions in S(tp) whose restriction to the [0, 1] 
form a basis for the space of all the functions in S(tp), restricted to the [0, 1]. 

When tp is the B-spline of order m (so N = m), then all polynomials of degree 
less or equal than m— 1 are in S(tp). Moreover, the set {1, x, x 2 , . . . , a;" 1-1 } is a local 
basis for S(tp), and the spectrum of T consists exactly of {1, 2 _1 , . . . , 2~( m ~ 1 )}. 

Now, if tp is not a B-spline then T could have some eigenvalue A different from a 
power of 1/2. If powers of 1/2 are associated to homogeneous polynomials, which 
functions in S(tp) are associated to an arbitrary eigenvalue A? Will the functions, 
associated to all the eigenvalues, provide also local bases of S(tp), or equivalcntly, 
will these functions reproduce the generator tpl 

In the case that A is a simple eigenvalue, Blu and Unser |BU02| and later Zhou 
Zho02' showed that A is associated to what they call a 2-scale A-homogeneous 
function, that is a function in the SIS that satisfies the relation h{x) — Xh(2x). 

However, to obtain a complete representation of the space it is necessary to con- 
sider the whole spectrum of T. This motivates the study of the spectral properties 
of T for a general refinable tp. This is achieved in this paper: we are able to com- 
pletely characterize the SIS in terms of functions associated to the spectrum of T. 
We prove that these functions provide a local basis of S(tp) (c.f. Theorem|3J). The 
advantage of this local basis is that it contains all possible homogeneous polynomi- 
als in the space, and those functions in the basis which are not polynomials, still 
preserve some kind of homogeneity. Furthermore this basis can be easily obtained 
from the spectrum of the finite matrix T. We are also able to prove that this matrix 
is necessarily invertible if the translates of tp are linearly independent. 

Definition 1. Let A S C, A ^ and r > 1 an integer. A function h is (2, A, r)- 

homogeneous if it satisfies the following equation: 

(1.2) £Q(-Ar fe M2- fc *)=0 a.e.. 

fc=0 ^ ' 

r is called the order of homogeneity, and A the degree. 
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If Ti C S(ip) is the span of all the (2, A, r)-homogeneous functions in S(ip), for 
any A 6 C and any positive integer r, we show that under the hypothesis of linear 
independence of the translates of ip, dim(7Y) = N+l, and that there is a basis of Ti, 
corresponding to the spectrum of T. More precisely, given a basis B = {v$, . . . , vn} 
of C N+1 that yields the Jordan form of T we associate to each vector v € B a unique 
(2, A, r)-homogeneous function in S(tp), where A and r satisfy v(T — AJ) r = 0. 

The first N of these functions are a local basis of functions in S(ip) restricted to 
[0, 1]. This allows to reconstruct the generator ip from the homogeneous functions, 
and gives a new representation for the functions in S(ip). 

Furthermore, we show that to each non-zero vector in the kernel of T, there 
corresponds a non-trivial linear combination of the integer translates of ip yielding 
the zero function. 



2. Notation 

Let if : M. — > C be a function supported in [0, N] satisfying 

We will often use an infinite column vector associated to ip, namely 

(2.1) = [..., <p(x - 1), ip(x), tp(x + 1), . . . ] . 

Let f(Z) be the space of all the sequences defined in Z. We say that the in- 
teger translates of ip are globally linearly independent, or linearly independent if 
SfcgZ a k l fi(' — k) = => ctk = Vfc, for any sequence a £ t(1). 
The subdivision operator associated to the mask Cfc is the operator 

(2.2) S c : £(Z) £(Z) defined by S c (a)j = ^ a % c 2l - v 

lG Z 

Note. The subdivision operator is sometimes defined in a different but equivalent 
way as S c (a)j — J2i<£Z aiC i~ 2i ' ^ ^ : ^(^) ~~ > ^(^) ^ s * ne °P era t or h( a )k = a -k, 
then hS c h — S c and therefore S c and S c share most of the properties. For a nice 
account of properties of the subdivision operator see jB.T02| . 

If L = L v is the double infinite matrix L = [c2i-j] i then the refinement 
equation can be written as 4>(x) — L<f>{2x). 

Using the matrix L, the subdivision operator l|2.2|l can be recast as: S c a = 
aL, a £ £{Z), where a on the right hand side of the equation is thought as an 
infinite row vector. Note that the scaling matrix T defined in Proposition^ is a 
finite submatrix of L. We will consider in our analysis the matrices M,Tq,T\, 
that are submatrices of T and are defined as: M — [c2i-j]i,j=i t ...,N-U Tq = 
[c2i-j]i,j=o,...,iV-ij Ti — [c2i-j]i ) j=i,...,Ar- That is, 



" c 










M 










cn 



Note that cq and cn must be non-zero, since Supp(y>) = [0, N). 

Now, if Y E ^(Z), define Y° and Y M as the restriction of Y to the indexes 
{0, . . . , N}, and {1, . . . ,N — 1}, respectively, i.e., 

Y° = (Y ,..., Y N ), Y M = (Yi, . . . , yjv-i). 
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Note. Throughout this paper, (L — XI) is considered as an operator on £(Z), defined 
by left-multiplication, (i.e. Y i — > Y(L — XI), where Y is a double infinite row 
vector). / is the identity operator acting on £(Z). By an abuse of notation, we will 
use the notation / for all identity operators, without distinguishing the space they 
are acting on. 

3. The point spectrum of L 

The following proposition, will show, how the spectral properties of L are related 
to those of T. The case r = 1 has been studied earlier by jCHMOOj . |JRZ98| . 
|Zho01llZno02| . 

Proposition 2. Let X G C. 

(1) Let Y G £(Z) andr G N, r > 1. IfY e Ker(L - A/) 1 *, then Y° G Ker(T - 

A/) r . Moreover, if X ^ 0, Y ^ and Y G Ker(i - XI) r , then Y° ^ 0. 
^ If v £ Ker(T — AJ) r and A ^ 0, £/ien £/iere exists an extension Y v G £(Z) o/ 
('i.e. Y® — v) such that Y v G Ker(Z, - A/) r . 

Proof. The matrix L (and therefore L — A/) can be decomposed in blocks as 



(3.1) L = 



R 





" 




'R-XI 








P 


T 


Q 


, (L - XI) = 


P 


T-XI 


Q 








s 










S-XI 



where we decompose Z as Z = A~ U A U A+, with = Zfl (-co, -1], A = 
Z n [0, AT] and A+ = Z D [N + 1, +oo), and 

i? = i| A - xA - P = i| A o xA - T = i| A o xA o Q = iU°xA+ <5 = iU+xA+- 

This block form of the matrix, is closed under multiplication. So if r > 1, r G N 
(3.2) 





' R r 












" (R-xiy 








u = 


Pr 




Qr 


, and (L-XI) r 




P \ 

r 


(T - A/) r 


Qr 










s r 












(s - xiy 


where P r = 




1 T k PR r - k - 1 and Q r = 


El 


Zl T k QS r ~ k 


-1 , and P, A 


and Q A are 



analogous, with the obvious changes. 



Note that the matrix S is upper triangular, with diagonal (0, 0, 0, ... ), and hence 
(S — XI) r is upper triangular, with diagonal ((— A) r , (— A) r , (— A r ), . . . ). 

Analogously, we observe that R is lower triangular with zeroes in the main 
diagonal, so (R — XI) r is lower triangular with diagonal ((— A) r , (— A) r , (~A) r , . . . ). 

If Y = (Y~, Y°, Y+), then Y(L- XI) r can be written as 

(3.3) (Y-(R - Xiy + Y a P*,Y°(T - XI) r ,Y Q Q* + Y+(S- XI) r ). 

So if Y G Kcr(L - A/)' r , then Y a G Kcr(T - A/) r . 

We now want to show that if Y G Ker(L - A/) r , A ^ 0, Y ^ 0, then Y° ^ 0. 

For this, let fc G Z be such that Yfe ^ 0. If < k < N, we are done. 
Assume that fco > N. Then, since Y(L — A/) r = 0, in particular, the fco element 
of this product is 0. But since A ^ 0, (S — XI) r is upper triangular with (— A) r 
in the diagonal, therefore the only nonzero elements of column fco of [L — XI) r are 
between and fco- Hence there has to be a fci, < fci < fco such that Y^ ^ 0. 
Again, if < fci < N we are done, otherwise we repeat the argument until kj is in 
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the desired interval. If fco < 0, the argument works in the same way, reversing the 
role of (S - A/) r and (R - AI) r . 

For the proof of partial assume that v G C N+1 , v G Ker(T — A/) r . We want to 
find an infinite vector Y G ^(Z), such that Y° — v and Y G Ker(L — A/) r . From 
equation 1(3.3(1 we know that if Y G £(Z), 

[Y(L-XI) r } + = Y Q Q* + Y+(S-\I) r , [Y{L-XI) r ]" = Y°P? + Y-(R-\I) r . 

Therefore, if Y G Kcr(L — \I) r , and Y° = v, then Y + and Y~ have to satisfy 

(3.4) Y + (S — XI) 7 ' = -vQ* and Y~(R- XI) r = -vP r x . 

Using again that (S — XI) r and (R — XI) r , are triangular, if A ^ 0, there are unique 
solutions for Y + and Y~ and they can be obtained recursively. □ 

The last proposition tells us that the elements of the spectrum of T are intimately 
related to those of L. But by the special form of T (see equation 1)2.3(1 ). we can 
actually use the (TV — 1) x (N — 1) matrix M to obtain the spectrum of T, as the 
following proposition shows: 

Proposition 3. Let A ^ £ C. 

(1) Letv = (v , ...,v N )e C N+1 andre N, r > 1. Then, ifv° G Ker(T-AI) r 
then v M = (vi, . . . , Ujv-i) G Ker(M — A/) r . Moreover, if X ^ Co, A ^ cjy, 
v° G Ker(T - XI) r , and v° ^ 0, then v M ^ 0. 

(2) if v AI — (vi, . . . , v N— i ) G Ker(M — A/) r and X ^ cq and X ^ cm, then there 
exists an extension v° G C^^ 1 of v , such that v° G Ker(T — XI) r . 

Proof. The proof is immediate by noting the special block-form of T. □ 



The case A = could not be handled with the methods of Proposition^ since the 
matrices R and S in (|3.1|) have zeros in the main diagonal. Instead, we need some 
results from the theory of difference equations which we present below |Hen62| . 

4.1. Difference Equations. Consider the linear difference equation with constant 
coefficients of order r 

(4.1) u y„ + uiy„ +1 H h u r y n +r = y = {y n } ne z, 

where Uk £ C,«o ^ 0,ti r ^ with characteristic polynomial P(x) = Y^k=o u k% ■ 

A solution to the equation ((4.1(1 is a sequence Y in £(Z), that satisfies 1(4.1(1 for 
all fceZ. A vector y — (yo, ■ ■ ■ , Vm) with m > r + 1 is a finite solution of 1(4.1(1 . if 
it satisfies ((4.111 for n = to n = m — r — 1. 

The space of solutions S C €(Z), has dimension r, and a basis of this space (the 
fundamental basis) can be written in the following way: 

Let ft, > 1 be an integer, d\, . . . ,dh arbitrary non-zero complex numbers with 
di dj if i =/= j. Let r 1; ...,r^ be positive integers. To each pair (e^,^), i — 
1, . . . , h we will associate a sequence — {a^} ke % defined as follows: Set r = 
7"i + • • • + 77j and ro = 0. Let < i < r — 1 and s = s(i),j — j(i) be the unique 
integers that satisfy ro + • • • + r s _i < i < r\ + ■ ■ ■ + r s , = i — J^t—o * rfe - 



4. The kernel of L 



(4.2) 



Define 




for \k\ > 

\k\<m- 
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where sg(fc) is the sign of fc. 

So, if P is the characteristic polynomial associated to 1)4.1(1 . consider the pairs 
{(di,r,) : where di is a root of P and r, its multiplicity}. The sequences {aife} fe6 2' 
i = 0, . . . ,r — 1 form a basis of S, the subspace of £(Z) of the solutions to (|4.1() . 

It is also known from the theory of difference equations, that every solution is 
determined unequivocally by any r consecutive elements of it. Hence, if y is a 
solution such that r consecutive elements are 0, then y is the zero solution. 

We will now associate to the pairs {(di,Ti) : i = l,...,h} the r x r matrix 

A = [<Hj]i,j=o r -i- Then (cf. Henrici Hcn62], pg. 214) 

h 

(4.3) det(A)= H (d, - d s ) r > +r ° Y[(n - 1)!!, 

l<i<S</l 1=1 

where 0!! = 1 and fc!! = kl(k - 1)! . . . 1!. Since d t 7^ d for i ^ j, det(A) ^ and A 
is invertiblc. 

Let us now consider a system of fc linear difference equations with constant 
coefficients of order r. 

(4.4) u i0 y n -\ h u ir y n+r = 0, i = l,...k, neZ, 

and let Pi be the characteristic polynomial of equation i, Pi(x) — Y)j— ugX 3 . 
Dchne T> = uf_ 1 {d : Pi(d) = 0} = {d±, . . . , d s }, and for each d G T> define r d = 
max{ri : Ti is the multiplicity of d in Pj}. Note that r d > 1 V<i S T>. We then have 
the pairs (di, r^) = (di, r^). Define the index of the system to be t — ^2 deT> r d < kr. 
Let I be the degree of the maximum common divisor p, of {P a , ■ ■ ■ , Pk}- Hence, 
P^x) = p(x)Pi{x), with degree Pj = r - i. (Note that t could be 0). With the 
above notation, we have the following proposition: 

Proposition 4. The space Sk of solutions to the system (14. 4|) has dimension £, 
where £ is the degree of the maximum common divisor of the characteristic polyno- 
mials. Moreover, if t is the index of the system 1)4.4(1 . and z is a vector of length t 
that satisfies 1(4.4(1 . then it can be extended to a sequence y z — {yj} e ^ solution of 
((4.4(1 and such that yj = Zj, j = 1, . . . , t. 

Proof. Let p be the maximum common divisor of Pi, . . . , P&, and let £ = deg(p). It 
is clear, that dim(5fc) > £. 

For the other inequality, consider the t x t matrix A = [a,ij]i j = o,...,t—i] with ay 
defined in ((4.2(1 for the pairs {(di,n)} defined above and t being the index of the 
system ((4.411 . Since di ^ dj for i ^ j, det(A) ^ by ((4.3(1 and A is invertible. 

Assume now that y G Sk, then y is a solution to all fc difference equations, hence 
there exist a 1 , . . . , a k vectors of length r, such that 

(4.5) A t a l = [mo, — , We— x]* 1 < » < fc, 

where A 1 is an t x r matrix whose columns are a fundamental system for equation 
i. Note that A 1 is a sub-matrix of A, whose columns correspond to some columns 
{ii, . . . ,i r } of A. 

Let now a % be vectors of length t, such that <S^ = whenever h $ ,i r } 
and a\ = a s , s = 1, . . . , r. Then we have for i, j = 1, . . . , fc 

(4.6) Atf = A ia * = [y , yt^f = A ja j = Aa J , 

and hence A(a l — — 0, for all i ^ j and therefore, by the invertibility of A, 
a 1 = oP , for all i ^ j. Therefore the only non-zero elements of a 2 : can be those 
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corresponding to the columns associated to the roots of p. Hence y is a linear 
combination of £ columns, and therefore dim(S'fc) < I. □ 

By noting that for the previous proof, we only used the first t coordinates of the 
infinite sequences, we have the following immediate Corollary. 

Corollary 1. If z is a vector of length t that satisfies l|4. 4JI . then it can be extended 
to a sequence y z = {yj}j e ^ solution of l|4.4|) and such that yj — Zj, j = 1, . . . ,t. 

4.2. The Ker(L). We can now return to our double infinite matrix L and look at 
the special case A = 0. As it turns out, the kernel of L is characterized by the 
vectors in the kernel of M. Since Co and cn are non-zero, the matrices T and M 
have kernels of the same dimension. Moreover, we have the following Proposition: 

Proposition 5. Consider the polynomials p e andp a of degree q — jv ~ 1 (we assume 
N to be odd) p e (x) = c + c 2 x + • • • + c 2q x q , p a (x) = C\ + c 3 x + • • • + c 2q+1 x g . 
Then dim(Ker(L)) = dim(Ker(M )) = degree(p), where p is the maximum common 
divisor of the polynomials p e and p . In particular, if dim(Ker(M)) > 0, p e and p D 
have a common root. Furthermore 

(1) For every Y £ Ker(L), Y ^ 0, we have Y M ^ and Y M £ Ker(M). 

(2) Conversely, for each v £ Ker(Af), v ^ 0, we have Y v ^ and Y v £ Ker(L). 

Proof. Let us observe first, that Y £ £(Z) is in the Kernel of L, if and only if Y 
satisfies the system of difference equations: 

f C V n + C 2 V n +l + h C 2q V n +q = 

\ cii>„ + c 3 v n+1 H h c 2q+1 v n+q = 0. 

Therefore, by Proposition^] Ker(L) is the subspace generated by the fundamental 
solutions associated to the roots of p, the maximum common divisor of p Q and p e . 
This shows that dim(Ker(L)) = degree(p). 

On the other side, if Y £ Ker(L), since (YL) M = Y M M we conclude that 
Y M £ Ker(M), and if Y M is the zero vector, then the solution Y of l|4~7jl has N- 1 
consecutive zeros, so Y — 0. Hence, if Y ^ 0, then Y M ^ 0, which proves (JIJ. 

To see that if v = (vi, . . . , ujv-i) satisfies vM — 0, then v can be extended, just 
note that the sequence v\, . . . ,vn-i must satisfy the difference equations system 
of order q — (we assumed N to be odd) given by (|4.7|) . Since the index t of 
the system |4~7I satisfies t < 2q = N — 1, and v is a non-trivial common solution of 
length N — 1, by Corollary 2] this solution can be extended in such a way that the 
extension satisfies both difference equations. This proves (J2J. 

From JIJ and J2J it is immediate that dim(Ker(L)) = dim(Ker(M)). □ 

Note. The fact that dim(Ker(M)) > implies that p e and p have a common 
root was proved under some minor technical conditions by Meyer |Mey91| . Related 
results can also be found in | JW93j . 

4.3. Invertibility of L. Propositions[21and|ni relate the spectral properties of the 
matrix M to the ones of the operator L. The next proposition shows a necessary 
condition for the independence of the integer translates of the function ip, in terms 
of the matrix M. 

Proposition 6. With the above notation, consider the following properties 
(1) {cp(- — k)} ke ^ are globally linearly independent. 
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(2) The operator L : — ► ^(Z), Y i — > YL is one-to-one, 

(3) The matrix M is invertible. 
Then and 0=^0). 

Proof. (Jl]=^|2l) Assume YL = 0. Define F(x) = Y<j)(x). Then we have F(x) = 

Y<j)(x) = YL<j){2x) = 0. Now, Y<f>(x) = => Y = 0, therefore Ker(L) = {0}. 

(HI EI) is a consequence of Proposition [5] □ 

Note: We do not know if either J2J or implies 

5. Homogeneous functions 

Assume now that Y <G Ker(i — XI) r , and define the function /i e <S(y>) as 
/i(x) = Y"^(a;). So /i is (2, A, r) homogeneous (c.f. I|l-2fl ). since it satisfies: 

= V(i - A/)VM = V [it Q (-A)*i r -*J *(*) 

= >' ft © j^r)) - £ (!) t-*)'^)- 

\fc=0 v 7 / fc=0 v 7 

We will denote by H(2, A, r), the space of all (2, A, r) homogeneous functions. 

Remark (1). As pointed out by the referee, a (2, A, r)-homogeneous function, is 
a particular case of a poly-scale refinable distribution. The concept of poly-scale 
refinable distribution is a generalization of refmabihty. See for instance [DD02 
ISun()5| . 

Remark (2). Note that if h G H(2, A, r) then h G H(2, A, s) for every s > r. There- 
fore the "order of homogeneity" will be defined by min{s : h G K(2, A, s)}. 

Remark (3). If ft is homogeneous (of any order) and A =/= 1, then ft(0) = 0. The val- 
ues of any homogeneous function of order r in (0, +oo), are completely determined 
by its values on any interval of the type [p^r, p), k G Z. (Analogously, the values 
on (—00, 0), are obtained from the values in any interval of the type (— p, — pW] ). 

Remark (4). In the case of order of homogeneity 1, (e.g. r = 1), h is a 2-scale 
homogeneous function as described in |Zho02| . 

Proposition 7. Assume {ip(- — k)} are linearly independent. Let <f> be as in (|2.1|l . 

If gi, . . . , g n 6 S(cp), gi = Y l <\>, then {<?i, . . . , g n } are linearly independent functions 
if and only if {Y , . . . , Y n } are linearly independent in £(Z). 

Proof. We observe that Y^=i ai 9 i = ai (Y l <ft) = a iY l ) 4 1 - This equa- 

tion, together with the linear independence of the translates of <p, tells us that 
OLigi = if and only if (J2 i o^F 1 ) = 0, which proves the desired result. □ 

Theorem 1. Assume that {ip(- — k)} are linearly independent. If h — Y(f>, (h G 
S(ip)), and h G 7i(2,A, r) then Vh = Y a G Kcr(T — A/) r . Reciprocally, if v E 
Ker(T — A/) r , then the function h = Y v cf> is in H(2, A,r). (Here Y v is the unique 
extension of v to a vector in Ker(L — A/) r by Prop. \Q) 

Proof. For the first claim, note that 
= E (i) (-^) r - k H2- k x) = J2 ([) (-^- k YL k Hx) =Y(L- XI) r <f>(x). 

fc=0 ^ 7 k=0 ^ 7 
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Then, Y (L — XI) r — and by Proposition v h G Ker(T - A/) r . 

For the converse first observe that if v = the result is trivial. Assume v ^ 
and v G Ker(T — A/) r , then by Proposition [fj] A ^ 0. Hence (by Prop.|2J) there is a 
unique extension Y v G Ker(L — A/) r , so h = is in 7i(2, A, r). □ 

5.1. Local basis of homogeneous functions. 

Theorem 2. Lei A be the set of eigenvalues of T , and let B = {vq, ■ ■ ■ , i>at} be a 
basis of C +1 that gives the Jordan form of T . Let Tt = (BasA ^ A c SQp)i where 
7i\(ip) = S(ip) : h G 7i(2,A, fc), for some k > 1},A G A. T/ien we have that 
dim(H) = AT + 1. 

Remark. Note that we can choose both vq = (1, 0, . . . , 0) and vn — (0, . . . , 0, 1) to 
be in the basis B, corresponding to the eigenvalues cq and cn respectively. 

Proof. If Vi G B, (0 < i < A), then v, G Ker(T - XI) k and ^ £ Ker(T - A/)* -1 , 
for some A G A, and fc > 1. So to each i?j G B, we can associate a unique pair 
(A, fc). Let us denote such Vi = v(X,k). (Note that by the previous observation, 
v = v(cq, 1) and v N = v(c N , 1)). 

After Theorem^] we can associate to each v(X, fc) a function h v i\^) m ^(2, A, fc)n 
S(ip). Furthermore, the functions {h v (x.k)}veBi are linearly independent. For this, 
observe that since the vectors in B are linearly independent, its extensions {Y v } are 
linearly independent in ^(Z), and therefore the functions {/WA,fc)}i>eB are linearly 
independent. 

One can see that if a finite number of functions are homogeneous for the same 
A, then a linear combination of them is also homogeneous for the same A. More 
precisely, X)"=o a ihi{X, fc;) — h(X, fc), where fc = max,-(fcj). 

Hence, if px, the characteristic polynomial of T, is factorized as: pt{x) = 
Y\\ea( x — ^) rA ; tncn dim(H\) = r\ and a basis of H\ is the set of (2,A,fc)- 
homogeneous functions associated to the vectors v G B, such that v = v(X,k), 
for some fc > 1. □ 

Theorem 3. Assume that tp satisfies 1)1.1(1 and has linearly independent integer 
translates. Let B = {vq, . . . ,vn} be as before, a Jordan basis for T, and let B be 
the (N + 1) x (A + 1) matrix that has the vectors Vi as rows. Let 

4>°(x) — [ip(x), ip(x + 1), . . . , tp(x + A)]* and h(x) — [ho(x), h\{x), . . . , h N {x)] t , 

where hi is the homogeneous function associated to the vector Wj. Then we have 

(i) h(x) = B<jP{x) x G [-1,1]. 

(ii) 4>°(x) = T(fP(2x) and h(x) = BTB~ 1 h(2x) x G [-1/2, 1/2] where BTB^ 1 
is in Jordan form. 

(iii) There exists a local basis ofS(ip) consisting of homogeneous functions. More- 
over, if ip has accuracy n, this basis can be chosen to contain the polynomials 
{!••'••••• ; !- 

Remark. Note that JuJ is a statement about the refinability of both cf>° and h, where 
the scaling matrix of h is the Jordan form of the scaling matrix of </>°. 

Proof. Since the support of ip is [0, N], ip(x + fc) = if fc ^ {0,1,..., A} for 
x G [-1,1]. Then 

hi{x) = Y l <P{x) = Vi$°(x) x G [-1, 1]. 
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So we have, 

(5.1) h(x) = B(j>°{x) x G [—1, 1]. 

This shows that for every interval / C [—1, 1] the functions {ip(x), p(x+l), ■ ■ ■ , tp(x+ 
N)} span the same space than the functions {ho, ■ ■ ■ , h^} when restricted to /. 

Since {ip(x),p(x + 1),- • • , p(x + N — 1)} is a local basis of S(ip), if vn has 
been chosen to be vn = (0, • • • , 0, 1), using equation (|5.1|l {h , ■ • • , ft-Tv-i} are a 
local basis for S(ip). Moreover, if ip has accuracy n < N, then we can choose 
vq, ■ ■ ■ , v n -i to be the eigenvectors associated to the eigenvalues {1, • • • , 2~' i+1 }, 
and hence {ho, ■ • • , htt—i} = {lj x,-- ■ , a;™ -1 , h n , ■ ■ ■ , /ijv-i}. 

This proves 0) and jmj. 

For using again that the support of <^ is [0,7V], it is easily seen that if 
x G [-i, |] then <£(x) = 7>(2x). Then we have B<£°(x) = BTB- 1 B^°(2x) x G 

hill- □ 

5.2. Generalizations. Throughout the paper "function" meant "measurable func- 
tion" . However, with the obvious modifications, all results hold for the case that ip 
is a generalized function or distribution. 

It is interesting to consider the generalization of the results to arbitrary dilations 
M > 1 in R, and in higher dimensions with arbitrary dilation matrices. It can be 
shown that most of the results are still true in these cases CHM05 . 

6. Examples for N = 3 

B-spline The simplest case of refinable functions are the B-splines. They are 
the (normalized) convolutions of the characteristic function of [0, 1] with itself. In 
particular the B-spline of order 3 is the refinable function that satisfies: 

(6.1) b(x) = h(2x) + h{2x - 1) + h(2x - 2) + hj{2x - 3). 

The B-splines are those functions, for which the accuracy is maximum and so 
coincides with the dimension of the matrix To, so in this case, the eigenvalues of 
To are 1 (for the constant functions), ^ (for the linear functions), and \ (for the 
quadratic functions). 

Daubechies D4 Daubechies wavelets, are those refinable functions of N coeffi- 
cients, that are orthogonal and provide the highest order of accuracy possible. 
(Note that the splines do not form an orthonormal basis). D4 satisfies: 

D 4 (x) = 1 + ^ D 4 (2x) + 3 + ^ D 4 (2x-1)+ 3 D A (2x-2)+ l ~ ^ D 4 (2x-3). 

D4 has accuracy 2 (it reproduces the constant and the linear functions) . In this case 
the matrix T has eigenvalues 1, \ and c = 1+ 4 ■ So a basis for span{Z? 4 (a;), D 4 (x— 
1), D 4 (x — 2)} x£ [ 0i i] is also given by span{l, x, h C0 (x)} xG [ ,i] where h co is the ho- 
mogeneous function associated to cq- 

(A, l)-Homogeneous functions are not enough In the two previous examples, 
we could always obtain a basis of span{/(a;), f(x — l),f{x — 2)} xe [ 0! i] just by using 
1-homogeneous functions. The following example is to illustrate, that even in the 
simple case of only 4 coefficients, it may be necessary to use homogeneous functions 
of order bigger than 1. Consider the function: 

(6.2) f(x) = i/(2a;) + \f&x - 1) + r 3 f(2x - 2) + ~f(2x - 3). 
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It can be shown that / has accuracy 1, and the eigenvalues of T are {1, |}. So in this 
case, span{f(x),f(x - 1), f(x - 2)} xe[0>1] = spaa{l,h {1/ 3 >1} (x),h {1/3>2 }{x)} xe[0>1] , 
where ^{1/3,1} is a 1-homogeneous function corresponding to the eigenvalue 1/3, 
and /i{i/3.2} is a 2-homogeneous function corresponding to the eigenvalue 1/3. 
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